KH 15D: a star eclipsed by a large-scale dusty vortex? 

P. Barge and M. Viton 

Laboratoire d'Astrophysique de Marseille, Traverse du Siphon B.P.8, 13376 Marseille cedex 
12, France; Pierre.Barge@oamp.fr, Maurice.Viton@oamp.fr 

ABSTRACT 

We propose that the large photometric variations of KH 15D are due to an 
eclipsing swarm of solid particles, trapped in a giant gaseous vortex rotating at 
~ 0.2 AU from the star. The efficiency of the capture-in-vortex mechanism easily 
explains the observed large optical depth. The weaker opacity at mid-eclipse is 
consistent with a size segregation and a slow concentration of the particles toward 
the center of the vortex. This dusty structure must extend over ~ 1/3 of an 
orbit to account for the long eclipse duration. The estimated size of the trapped 
particles is found to range from 1 cm to 10 cm, consistent with the grey extinction 
of the star. The observations of KH 15D support the idea that giant vortices can 
grow in circumstellar disks and play a central role in planet formation. 

Subject headings: circumstellar matter - hydrodynamics - planetary systems - 
stars: pre-main-sequence 

1. Introduction 

KH 15D is a pre-main sequence (PMS) K7 star located in the young star cluster, 
NGC 2264. Its large photometric variations are periodic and possibly interpreted in terms of 
eclipses by an opaque feature orbiting inside a circumstellar disk seen nearly edge-on (Herbst 
et al. 2002). It is likely connected to the UXOri star family, characterized by large photo- 
metric and polarimetric variabilities (eg. Natta et al. 1999), but sets apart from them with 
two important differences: no infrared excess in the disk and no significant color changes 
during the eclipses, both indicating particle sizes much larger than optical wavelengths. As 
the VRI colours of the star are consistent with its spectral type (Hamilton et al. 2001), the 
selective absorption is A v < 0.2, a value indeed in agreement with the foreground interstellar 
reddening, E B _ Y ~ 0.06 - 0.07, of stars in NGC 2264 (Perez et al. 1987; Sung et al. 1997). 
On the other hand the grey absorption, A , during the out-of-eclipse observations can be 
estimated indirectly using a well known relation (Barnes et al. 1978) between the visual 
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brightness parameter and the intrinsic colour (V — R)q. Using the absolute visual magnitude 
deduced from the distance modulus of the cluster, this relation provides a value of the stellar 
radius consistent with that derived by Hamilton et al. (2001) and with the PMS status only 
if the grey absorption is ~ 0.6 — 0.7 magnitude. 

The five eclipses observed since the discovery of KH15D (Kearns and Herbst 1998) 
provide us with the main characteristics of the occultation: a relative depth ~ 0.95, an 
approximate duration of 18 days and a period of 48.36 days. The absorbing cloud is deduced 
to orbit at ~ 0.2 AU from the star and to extend over 1/3 of an orbit (Hamilton et al. 2001). 
Such an occultation completely differs from a standard planetary transit and also displays 
other important peculiarities: (i) ingress and egress are remarkably steep, indicating an 
occulting clump with a very sharp edge; (ii) the minimum is not a flat land but contains a 
central bump with small scale features and a reverse peak close to mid-eclipse; (iii) when 
comparing two successive occultations, the reverse peak seems to shift from one side of the 
mid-eclipse to the other; (iv) since the first series of observations, the occultation has clearly 
widened and deepened in secular fashion. 

The various models suggested by Herbst et al. (2002) to explain these observations 
involve a companion, planet or brown dwarf, able to "shepherd" gravitationally an extended 
swarm of solid particles (warped density wave or ring arc) which is responsible for the eclipses. 
Possibly, these models could account for a number of the observed characteristics but, so 
far, none of them has been developped far enough to reproduce the observed light-curve and 
no observation has confirmed the presence of a "shepherd" . 

In this letter, we propose a new model in which KH 15D is periodically occulted by a 
dusty anticyclonic vortex persisting for a long time in a gas disk surrounding the star. The 
particles are confined inside the vortex by the gas friction and no gravitational assistance is 
needed to keep them together into a clump. The swarm of the trapped particles is organized 
by the gas dynamics with an absorption profile that pretty well explains most properties 
of the observed light-curve. The gaseous vortex associated to this dust swarm covers a 
wide azimuthal extent (see Fig.l) which seems quite puzzling to explain but looks like those 
produced in recent numerical simulations (Klahr and Bodenheimer 2003). Our model is 
based on a scenario proposed by Barge and Sommeria (1994,1995; thereafter BS94,BS95) in 
which the formation of planets begins inside persistent gaseous vortices; it relies on a number 
of recent works on the formation and evolution of vortices in circumstellar disks (Li et al. 
2001; Klahr and Bodenheimer 2003). 



- 3- 



2. The proposed model 

Our model starts from a circumstellar disk of gas and solid particles in which the mech- 
anisms of planetary formation are still at work. The gas is assumed to spread into a flared 
disk following a simple hydrostatic equilibrium in the vertical direction with a scale height 
H — Cs/Q, where C$ is the sound speed and Q is the Keplerian orbital frequency. Protoplan- 
etary disks are likely turbulent during a span of their live and could host large scale and long- 
lived vortices growing from turbulence by analogy with what happens in two-dimensional 
fluid dynamics, in which organized structures are known to emerge from random turbulence 
in rotating shear flows. Vortices spinning like the shear flow are robust and merge one an- 
other while those with the opposite sign are laminated by the shear (BS95 and references 
therein). However, such vortices may also result from some specific instabilities as for exam- 
ple: (i) the Rossby Wave instability investigated by Li et al. (2001) which requires a strong 
local maximum in entropy and density; (ii) a global baroclinic instability described recently 
by Klahr and Bodenheimer (2003) and which arises from a natural radial stratification of 
the gas flow. 

A number of numerical simulations show that large scale vortices can grow and survive 
for many rotation periods and look like either elongated structures stretched by compressibil- 
ity effects (Godon and Livio 2000) or a single vortex dominating the whole disk (Klahr and 
Bodenheimer 2003). Highly elongated vortices are found also in the inner regions of MHD 
accretion disks (Tagger and Pellat 1998, and private communication) or as exact solutions 
of the incompressible Euler equation with different aspect ratios (Chavanis 2000). On the 
other hand, the radial extent of vortices cannot exceed the thickness of the disk (R < H) 
since the velocity of the vortex, RQ, must be less than the sound speed C$ = HQ (BS95). 
Larger vortices would be destroyed by energy losses due to sound waves and density waves. 
The solid particles embedded in the gas disk are submitted to a friction drag depending on 
the mean-free-path A of the gas molecules relative to the particle size (or radius) s. In the 
inner disk region (r < 0.2AU) and for particles larger than a millimeter, the mean- free-path 
is less than size and the drag, caused by particle wake, reaches the Stokes regime (in the 
limit of particle Reynolds numbers Re p < 1) with a stopping-time (Chavanis 2000): 

Ts = 8s2p0 , (1) 
9o"g as f2A 

where a gas is the gas surface density and p — 2g/cm 3 is the density of material particles 
are made of. In contrast, the Epstein regime is reached for smaller particle size s < 9A/4. 
The particles are only submitted to star attraction and gas drag; their dynamical evolution 
depends on the single non-dimensioned friction parameter r s = QT S : (i) the lightest particles 
(t s <C 1) come at rest rapidly with the gas and are driven by the flow; (ii) the heaviest 
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particles (r s ^> 1) are nearly unaffected by the gas motion and keep a quasi-Keplerian 
motion. 

The capture of particles by a gaseous vortex has been explored first by BS94 and BS95 
using a simple model in which the velocity field is made up of concentric epicycles inside 
the vortex, while it matches a Keplerian flow at large distances. They found that such 
an anticyclonic vortex can capture and concentrate dust particles very efficiently: (i) light 
particles penetrate into the vortex and stop on streamlines close to the edge where they 
slowly shift toward the core; (iii) optimal particles, with r s ~ 1, sink deeply into the vortex 
and reach core streamlines. These results were confirmed for other velocity fields (Chavanis 
2000; De la Fuente et al. 2002). The vortex induces a segregation of the trapped particles 
(or a size sorting if the particles have the same composition) following the value of their 
friction parameter; this is illustrated in a number of numerical simulations (De la Fuente et 
al. 2002). 

The ultimate reason for this dynamical behaviour lies in the sign of the vortex rotation. 
Indeed, in a reference frame rotating with the vortex center, the Coriolis force can overcome 
the centrifugal force and pushes the particles toward the core if the vortex is an anticyclone, 
whereas both forces are conspiring to eject the particles for a cyclone. This capture-in- 
vortex mechanism is a very efficient one and results in strong density enhancements inside 
the vortex, by at least two order of magnitude in ~ 200 rotation periods. The capture rate is 
estimated under the assumption that the particles are continuously renewed near the vortex 
orbit due to the inward drift under the systematic headwind drag (Weidenschilling 1977, 
BS95). One obvious and important consequence of this density enhancement is that, inside 
a vortex, particle growth is made easier and will couple with confinement and segregation. 
On the other hand, dust is depleted from the region inside the vortex orbit as the particles 
either are feeding the vortex or are falling to the star under the systematic drift. 

Inside the vortex, the trapped particles are also submitted to a background small scale 
turbulence which makes them diffuse and tends to reduce their global concentration. Chava- 
nis (2000) investigated this question in terms of a diffusion equation in an idealized circular 
vortex and derived a time dependant solution for the surface density inside the vortex: 



°* « jg^W) exp 



(r — kr ) 2 



(2) 



the initial state being a delta function centered at r ; l d is the diffusion scale length and 
k = exp(— t/T capt ) in which T capt is the characteristic time for a particle to reach the center 
of the vortex. In the case of light particles and very elongated vortices T capt ~ 2q/(3Qr s ) 
and Id ~ ^/a v /qT s R v where R v is the vortex radius and q its aspect ratio; a v is the non- 
dimensioned parameter measuring the small scale turbulence efficiency inside the vortex. 
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In order to estimate the size of the trapped particles, we will choose a standard model 
of nebula, the minimum mass solar nebula, in which the surface densities (both for gas and 
particles) and the temperature are the decreasing power laws r -3 / 2 and r^ 1 / 2 , respectively; 
at 1 AU the densities are set to 1700 gem -2 for the gas and 20 gcm~ 2 for the particles, 
whereas the temperature is assumed to be 280 K. At 0.2 AU from the star and with our 
numerical values, the optimal size for particle capture is s op t — cm, 
i.e. pebble size. Lighter non optimal particles with r s < 1 can be also captured by the 
vortex but remain trapped at the periphery of the vortex. These particles remain in the 
Stokes regime, at 0.2 AU from the star, as far as their size is larger than the critical size 
s c = 9A/4 ~ 0.3 mm. 



3. Fitting the light-curve 

We assume the disk is seen nearly edge on, under an inclination i less than the flaring 
angle of the gas disk. % must be small enough for the line of sight to cross the vortex in 
its vertical extent, but also large enough to avoid the prohibitive optical depths near the 
mid-plane layers. Of course, this inclination has to be consistent with the out-of-eclipse grey 
extinction A ~ 0.7 magnitude deduced in Sect.l. In its motion around the star, the vortex 
periodically crosses the line of sight and the optical depth r is a function of time varying 
from A (out of eclipse) to A + 3.2 during the eclipses. Assuming the absorption is due to 
spherical particles with sizes larger than a millimeter, the opacity k reduces to ns 2 n and the 
optical depth r = J ndl is the familiar dust column density, where n is the local number 
density along the line of sight. 

During an occultation r is found to increase inside the vortex, from center to edge, then 
to fall down very steeply to reach the out-of-vortex level. This "opacity-curve" has a peculiar 
shape that our model can easily reproduce. 

First, one can guess that a circular blob of matter in which the opacity «(r) is radially 
increasing can conveniently mimic the optical depth profile. Then, we noticed that a simple 
linear dependance, k oc 1 + a r, with a cut off at the outer boundary permits to fit pretty well 
the data (Fig. 2). In our fit the optical depth has been computed following two steps: (i) an 
integration of n along a path crossing an equivalent circular vortex (ECV) of radius H and 
located at a distance r or b from the star; (ii) a circular anamorphosis by a factor q ~ 17 (the 
result of an integration along a path at angle u is reported to the direction of observation 
at angle q x u). Figure 1 illustrates the ECV model with a sketch of the integration path 
and the density contours of the vortex based on the linear opacity law and the appropriate 
anamorphosis. 
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This absorption profile corresponds to a shell-like distribution of the solid particles 
inside the vortex, reminiscent of the spherical structures observed in Hll regions. In our 
model such a density profile originates in the dynamics and the segregation of particles 
trapped in a vortex as described above. 

Surprisingly, a simple cut off at the vortex outer boundary can easily mimic the very 
steep fall during ingress or egress, consistent with an occultation by a sharp edge as proposed 
by Herbst et al. (2002). In fact, a cut off cannot receive true physical justification and the 
resulting profile cannot connect smoothly the out-of-eclipse level. A more realistic fit of 
ingress and egress directly arises from our model. Indeed, light sub-optimal particles are 
preferentially captured in the outer part of the vortex where they accumulate and slowly 
diffuse toward the core. The fit is realized thanks to the Gaussian density profile presented 
in Eq. 2 with the following values of the parameters: the center of the Gaussian is close to 
the vortex boundary (at a distance from the core ~ R v = q H); the 1/e half-width of the 
Gaussian is given by / dv /(l - ¥) ~ 0.07i? v . This is possible using the following assumptions: 
(1) the particles are renewed by capture every ~ 10 rotations of the vortex and have a friction 
parameter r s ~ 0.01, consistent with a capture at the vortex periphery and a slow drift toward 
the core in a time scale T capt — 2q/ (3f2r s ); (2) the parameter measuring the turbulence inside 
the vortex is a v ~ 0.01. As a result, the correponding size of particles confined in the outer 
regions is of the order of a centimeter. 

4. Discussion and conclusion 

The proposed model can account for the main characteristics of KH 15D's observations. 
It is able to reproduce the light-curve and permits to fit pretty well the data. 

(1) The period and duration of the eclipses may be reproduced by the rotation around 
the star of a giant swarm of solid particles, trapped in a persistent gaseous vortex. 

(2) During totality, the rise up of the flux can result from a lower opacity of the central 
regions, possibly coming from a particle size segregation (the larger the particle, the deeper 
and the faster the confinement inside the vortex). 

(3) The estimated size of the trapped particles ranges from 1 to 10 centimeters, consistent 
with the grey absorption observed in KH 15D. The presence and growth of bigger particles 
is very likely due to the high densities and the low relative velocities. This is in agreement 
with the expected predominance of large grains in PMS disks which seems also required to 
explain sub- millimeter and millimeter disk observations (eg. Natta et al. 1999). 
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(4) The shift of the reverse peak observed in successive transits could be explained by 
gas rotation inside the vortex, the inner zone rotating at half the orbital frequency similarly 
to what happens in smaller vortices (Chavanis 2000). 

We can also speculate that: (i) "secular" changes of the light-curve could result from 
a dynamical evolution of the swarm of particles; (ii) the wavy behaviour at egress could 
correspond to a turbulent wake of the vortex. 

One of the main assumption in our model is that the disk is seen nearly edge-on, at an 
angle % ~ 2 — 3°, and with an optical depth A ~ 0.6 — 0.7 magnitude during the out-of-eclipse 
periods. This is a possible situation if the dust particles have grown and settled down to the 
mid-plane forming a flat sub-disk dominated by large grains. The predominance of large 
particles in KH 15D is, indeed, completely consistent with our results. However, we want 
to stress that no numerical simulation is, nowadays, able to reproduce the huge azimutal 
extent of the vortex required to fit the observations. Further hydrodynamical simulations 
and theoretical modeling of the particle dynamics are necessary to confirm and improve our 
model. 

Finally, KH 15D appears a good target to test planetary formation models. This is 
presently the case with the scenario in which vortices can persist in a gas disk acting as 
traps for dust (protected against losses into the star) and as wombs for planetesimals. Of 
course, further observations of this object are required to constrain more firmly the model. 
The possibility to observe large scale structures in circumstellar disks has been recently 
investigated by Wolf and Klahr (2002), extrapolating the capabilities of the future giant 
interferometer ALMA. KH 15D and other similar objects are good opportunities for the 
future space missions like CoRoT, Kepler and Eddington. 

A further piece of information was recently provided by Hamilton et al. (2003) using 
high resolution UVES spectra which confirmed that KH15D is a weak-lined T Tauri star 
surrounded by an accretion disk with, possibly, a collimated bipolar jet. Such configuration 
could favor the accretion-ejection model of Tagger and Pellat (1998) in which persistent 
vortices can form in the inner region of magnetized disks. 

While submitting this letter we learned that H. Klahr also claims the existence of a 
vortex in the disk of KH15D to explain the peculiar light-curve of this object. His work 
seems complementary to ours since he is focusing on the way a giant vortex forms and 
persists in hydrodynamical simulations, and not on the trapping of solid particles in gaseous 
structures. 

We wish to thank W. Herbst and CM. Hamilton for communicating us the KH 15D data 
plotted in Fig. 2. We are also grateful to A. Llebaria for his help in handling the data. 
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Fig. 1. — The banana-like dusty vortex used to model the eclipses of KH 15D. Right: Sketch 
of the equivalent circular vortex (ECV) of radius H centered at a distance r orb from the 
star. The optical depth obtained by integrating the opacity along a path crossing the ECV 
at angle u is reported to the direction of observation at angle q x u. Left: Approximate 
density contours (0.15 to 0.9 in steps of 0.15) corresponding to the linear opacity law (see 
text), a crude model surprisingly representative of the observed extinction curve. Note that 
the ratio H/r ov \, = 0.2 used in this figure is much larger than actual. 
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Fig. 2. — Eclipses of KH 15D in J c -band magnitudes, phased with the period of 48.36 days. 
The plotted points are weighted means per observation night from the 2001—2002 observation 
campaign (W. Herbst and C. Hamilton, private communication). Dashed line is a fit obtained 
with the linear approximation k oc 1 + a r and a ~ 10. Solid line is the best fit obtained 
assuming k oc b r (b ~ 5) if r < H, and k given by the Gaussian model described in the text 
if r > H. 



